Detecting entangled states in graphene via crossed Andreev reflection 
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Shot noise cross-correlations across single layer graphene structures are calculated with insulators 
separating a superconducting region. A new feature of specular crossed Andreev reflection comes 
into play due to the unique band structure of graphene. This gives rise to a rich structure in the states 
of the electric current flowing across the graphene sheet. We identified a parametric regime where 
positive shot noise cross-correlations of the current appear signifying entanglement. In contrast to 
previous proposals the sign of the cross-correlations can be easily tuned by the application of a gate 
voltage. 

PACS numbers: 72.70.+m, 73.23.-b, 74.45. +c, 03.65.Ud 



I. INTRODUCTION 



tion applications^. 
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Shot noise is defined as the temporal fluctuation of 
electric current in a non-equilibrium set-up^. When the 
shot noise cross-correlations between two regions turn 
positive it signals the presence of electronic entangled 
states^. In order to make use of these correlations for 
quantum information purposes one would need to spa- 
tially separate the electrons without destroying the en- 
tanglement^. This is ideally detected by entangled elec- 
trons traversing different wires^. The quantum correla- 
tions can be provided by Cooper pairs present in super- 
conductors, which is the most entangled state found in 
nature. 

To intuitively understand how shot noise contains the 
signature of entanglement we resort to statistics. Shot 
noise cross-correlations are defined as cross-correlations 
of current fluctuations across two distinct regions. Ab- 
sence of entanglement leads to positive correlations for 
photons (bunching) and negative for electrons (anti- 
bunching) . The observation of positive shot noise correla- 
tions for electrons is a signature that they are in an entan- 
gled state. This has been most famously predicted in nor- 
mal metal-superconductor-normal metal structures'^, 
but it has not yet been experimentally demonstrated. 
An earlier experimental attempt in a two dimensional 
electron gas beam splitter connected to a superconductor 
could not arrive at any definite conclusion possibly due to 
the low tunability of these devices. In this work we inves- 
tigate what happens to the noise cross-correlations when 
a single layer of graphene replaces the normal metal or 
semiconductor. Our motivation comes from the follow- 
ing fact. In contrast to a normal metal, the energy of 
transported electrons can be very efficiently controlled 
in a graphene layer via the application of a gate voltage 
thus being much more amenable to experiments. This 
was demonstrated in Ref£ where it was shown that the 
Josephson current could be very efficiently tuned via the 
application of a small gate voltage. We expect that this 
characteristic of graphene structures will facilitate the 
observation of entanglement in solid-state systems, thus, 
opening the way for their wider use in quantum informa- 
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FIG. 1: An overview of the setting from the top. Two in- 
sulating layers of graphene (Gi's) on either side of the su- 
perconducting graphene layer (Gs). Voltages Vi and V2 are 
applied to either ends of the normal graphene layers (Gn's). 
Schematic of specular crossed Andreev reflection is also de- 
picted. Incident electron at angle 9 (IE). Reflected electron 
at angle — 9 (RE). Andreev reflected hole at angle 9a (AR). 
Specular Andreev reflected hole at angle —9a (SAR). Elec- 
tron like quasiparticle (ELQ). Hole like quasiparticle (HLQ). 
Crossed Andreev reflection at angle 9a (CAR). Specular 
crossed Andreev reflected hole at angle —9a (SCAR). Elec- 
tron co-tunnelling at angle 9 (EC). 



II. MODEL 

Graphene is a monatomic layer of graphite with 
a honeycomb lattice structured that can be split 
into two triangular sublattices A and B. The elec- 
tronic properties of graphene are effectively described 
by the Dirac equation^. The presence of isolated 
Fermi points, K + and K_, in its spectrum, gives 
rise to two distinctive valleys. In this work we 
deal with a normal-insulator-superconductor-insulator- 
normal (NISIN) graphene junction. We consider a sheet 
of graphene on the x-y plane. Superconductivity is in- 
duced via the proximity effect, where a normal supercon- 
ductor at close range on top of the sheet generates the 
desired superconducting correlations. In Fig. Q] we sketch 
our proposed system. The superconducting region is lo- 
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FIG. 2: Energy-Momentum diagram to explain specular 
crossed Andreev reflection where the regions of normal (N) 
graphene and superconducting (S) graphene are as indicated, 
(a) Ef > A regime where Andreev and crossed Andreev re- 
flection occur in the same band, and (b) Ef -C A regime 
where Andreev and crossed Andreev reflection occur in a spec- 
ular fashion across the bands. If, E > Ef > A, where E is 
the energy of the particle measured from the Fermi level, then 
also specular Andreev reflection occurs. In all calculations we 
are in the regime, where E < A. 



cated between < x < L, while the insulators are located 
on its left, — d < x < 0, and on its right, L < x < L + d. 
The normal graphene planes are to the left-end, x < —d, 
and to the right-end, x > L + d. 

As we will demonstrate in the following there are 
additional processes occurring at the normal graphenc- 
superconducting graphene-normal graphene junctions 
than those seen at normal metal-superconductor-normal 
metal junctions^. These are local specular Andreev re- 
flection and crossed (non-local) specular Andreev reflec- 
tion. Importantly, Andreev reflection in graphene can 
switch the valley bands, i.e., conduction to valley, see 
Fig. [21 This process is known as specular Andreev re- 
flection^, explained in Fig. [21 In the process of normal 
Andreev reflection, an incident electron from the normal 
metal side is reflected as a hole which retraces the tra- 
jectory of the electron. In specular Andreev reflection, 
the reflected hole follows the trajectory which a normally 
reflected electron would have. In this work we see, in ad- 
dition to this, the possibility of specular crossed Andreev 
reflection, where a hole is reflected at the other lead but 
in a specular fashion (see Fig. [I]). 

For a quantitative analysis we describe our system by 
the Dirac Bogoliubov-de Gennes equation that assumes 
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where E is the excitation energy, A is the superconduct- 
ing gap of a s-wave superconductor, \& is the wavefunction 
and : represents 4x4 matrices. In the above equation 
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H± = -ihv F (a x d x ±(jydy) + U (2) 



Here %, vf (set equal to unity hence forth) are the 
Planck's constant and the energy independent Fermi ve- 
locity for graphene, while the a's denote Pauli matrices 
that operate on the sublattices A oi B. U is the elec- 
trostatic potential which can be adjusted independently 
via a gate voltage or doping. We assume U = 0, in the 
normal regions, while U = Vi, i = 1, 2, in cither insulat- 
ing regions and U = —Uq in the superconductor. The 
subscripts of Hamiltonian ± refer to the valleys of K + 
and K- in the Brillouin zone. T = —r y <g> a y C, (C being 
complex conjugation) is the time reversal operator, with 
t being Pauli matrices that operate on the ± space and 
/ is the identity matrix. 

Let us consider an incident electron from the normal 
side of the junction (x < —d) with energy E. For a right 
moving electron with an incident angle 9 the eigenvector 
and corresponding momentum reads 



4>l = [1, e % \ 0, 0} 1 e lp cos6te , p e = (E + E F ). (3) 

A left moving electron is described by the substitution 
8 — ► 7T— 9. If Andrccv-rcflcction takes place, a left moving 
hole is generated with energy E, angle of reflection 9a 
and its corresponding wave-function is given by 

i\) h _ = [0,0,1, e-^Apg-y-cos^* p h = (E-E F ). (4) 

The superscript e (h) denotes an clcctron-likc (hole- 
like) excitation. Since translational invariance in the 
y-direction holds the corresponding component of mo- 
mentum is conserved. This condition allows for the de- 
termination of the Andreev reflection angle 9a through 
p h sm(0 A ) = p e sin(0). There is no Andreev reflec- 
tion and consequently no sub-gap conductance for an- 
gles of incidence above the critical angle 9 C = sin~ 1 (|i? — 
E F \/{E + Ep))- In the insulators, — d < x < and 
L < x < L + d, the eigenvector and momentum of a right 
moving electron are given by 

f|T ipj cos 9 io x 







with i 



= 1,2 while a left moving hole is described by 
[0, 0, 1, e- ie *} T e~ ip " cose ^ jP ^ = {E-E F 



Vi), 
(5) 



Vi). 
(6) 

On the superconducting side of the system, (0 < x < 
L), the possible wavefunctions for transmission of a right- 
moving quasiparticle with excitation energy E > read 
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with q e = (E F + Uq + ^E 2 - A 2 ) and q h = (E F + 
Uq — v 1 E 2 — A 2 ). In the sub-gap regime the quasi- 
particle wave-vectors have a small imaginary component 
as q e / h = E F + Uq ± l/£, where £ = 1/A is the coher- 
ence length. The coherence factors are given by u = 

yJ(i + y/T=&jw)/ 2 , v = ^(l-yr^AW)^. 

We have also defined 8 + = 9$, 9~ = n — 9$. The trans- 
mission angles 9$ for the electron-like and hole- like quasi- 
particles are given by q a sin 9g = p e sin 9, a = e, h. In the 
following we limit ourselves to the regime where Uo 7S> A, 
such that the mean field conditions for superconductivity 
arc satisfied. The trajectory of the quasi-particles in the 
insulating region arc defined by the angles 9io and 9^. 
These angles are related to the injection angles by 

sin 6 i0 / sin 9 = (E + E F )/(E + E F - Vj), 
sm9f /sm9= (E + E F )/(E - E F + V t ). (8) 

Here, we adopt the thin barrier limit defined as, 
9io,9^ and d — > 0, while Vi — > oo, such that p^d^p^d — > 
Xi- To solve the scattering problem, we match the 
wavefunctions at four interfaces: ip\x=-d — ^ l ii\x=-d 1 

lpll\x=0 = *s|x=0, *s|x=i = lp2l\x=L, and 4>2l\x=L+d = 

ip\ x =L+d, where, starting with normal graphene at left, 
V = V4 + s\\r- + stfyt i/> a = Pi^ e u+ + qii>tj_ + 
miipii+ + n^j_,i = 1,2, <J S = + q s V e s _ + 

m s^s+ + n s*s_, and finally for normal graphene at 
the right, ip = s^|-0+ + Si2"0+- Solving these equations 
leads to the amplitude of Andreev reflection sf[, normal 
reflection s\\, amplitude of electron co-tunnelling (EC) 
sf|, and of crossed Andreev reflection (CAR) sfrj. 
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FIG. 3: (a) Crossed Andreev reflection and (b) electronic co- 
tunneling as function of the superconducting length, L/£ . In 
both figures, xi = ~X2 = tt/ 4 > Uo = 1000A and E = 0.15A. 



III. RESULTS 
A. Specular crossed Andreev reflection 

The first issue we tackle is the non-local conductance. 
Similar caluclations, but for bipolar structures, were per- 
formed in Ref . [l4| ■ It is defined as the conductance in the 
right lead when both superconduction region and right 
graphene layer are grounded, while a voltage is applied 
to the left graphene sheet. The non-local conductance is 
given by the difference between the crossed Andreev and 
electronic co-tunneling currents in the absence of a bias 
at right, where G = Gcar — GeCi within 

G C ar = I" d9 cos 9 A \sf 2 \ 2 , G EC = /"^cos^l 2 . 

. (9) 

In the following figures all the quantities arc in their di- 
mensionless form with the superconducting gap set to 
A = I. The other energy parameters are expressed in 
terms of A. In Fig. [3] we plot the non-local CAR and 
EC current as function of the length of the supercon- 
ducting region. We differentiate between two regimes. 



First, for E F ^> A there is absence of interband non- 
local electron-hole transmission, denoted as the crossed 
Andreev regime (Fig. [2|). Second, for E F <C A non- 
local interband electron-hole transmission is permitted 
giving rise to the specular crossed Andreev regime. As 
function of the length we find that both non-local coef- 
ficients vanish for large values. However, while the EC 
current decreases almost monotonically from a peak at 
L <C £ to vanishing for L ^ £, the CAR current is maxi- 
mum around L ~ £, and it vanishes for the extreme cases 
L <C £, L S> £. In Fig. [5J we plot the crossed Andreev 
current for normal transmission (left) as well as specular 
reflection (right). We observe that the specular CAR cur- 
rent might dominate the normal current in the E <C A 
regime. One very interesting fact, which is partly seen 
in NS graphene junctions, is that, just like the specular 
Andreev reflection seen there, here too the crossed specu- 
lar Andreev reflection is reduced to vanishing at E ~ A, 
but the normal crossed Andreev current is marginally re- 
duced at E ~ A. However, the non-local conductance 
(see Fig. [5]), is dominated by electron co-tunnclling. It is 
also periodic as function of the strength of the insulating 
barrier's \i' s ( n °t plotted here)^. 
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FIG. 4: (a) Crossed Andreev reflection and (b) specular 
crossed Andreev reflection. Non-local conductance for a 
NISIN graphene based structure as function of the electronic 
energy for (c) normal and (d) specular reflection cases. In all 
figures, xi = — X2 = tt/4, and Uo = 1000A. 



B. Shot noise cross-correlations 

Next we calculate the shot noise cross-correlations, 
which is the main focus of our work. For that we first 
have to derive an expression for the shot noise in multi- 
terminal settings^ applied to graphene. The fluctuations 
of the current away from the average is termed noise. A 
general expression for current fluctuations between any 
two arbitrary leads is given by 

N tj (r) = (AI t (t)Al) (t + r) + Al) (t + r)Al) (t)) , (10) 

where AIi(t) = Ii(t) — (Ii(t)). The Fourier transform of 
Eq. |10]) gives 



N ij {w)d(w + w') = (Aliiw^AIjiw) + i j (w)AI l (w')). 

(11) 

For simplicity we consider the experimentally feasible 
zero frequency noise limit, where displacement currents 
are absent. The current operator is given by 



Ii(w = 0) 
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FIG. 5: Noise cross-correlations as function of the gate voltage 
(x = Xi) applied to the left insulator. The right insulator is 
fixed at gate voltage \2 = 0, while (To = 1000A and V2 = 
0.2A. 



with A kl ,u{i,a) = Sik5ii8 ai 8 a s - s" fc 7t sf/, where Greek 
indices denote the nature (e for electrons, h for holes) 
of the incoming/outgoing particles with their associated 
charges q a , while Latin indices I, k identify the graphene 
sheets and au denotes annihilation operator for a parti- 
cle in lead I with charge 8. From Eqs. (fTT|) and (fl~2|l the 
zero frequency noise cross-correlations between the cur- 
rents at left and right normal graphene sheets (Gni, Gnz) 
become^ 



12 : 



E 



k, I 6 Gn 1 ,Gn 2 ,Gs 
a, /3, 7, S £ e, h 



dOcosO J dEAk-y,ts(l,a) 
xA lSM (2,/3)f kj (l-f ls ) (13) 



fk-y is a Fermi function for particles of type 7 in graphene 
sheet k. 

In the limit L <C £ Andreev and cross-Andreev re- 
flection vanish, which implies that in this limit noise- 
correlations are negative^. In the limit L 3> £ both 
non-local currents vanish leading to vanishing noise cross- 
correlations. However, it is the length in-between these 
limits where noise not only becomes substantial but also 
can change sign. In Fig. [5] we plot the shot noise cross- 
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correlations as function of the gate voltage, which tunes 
the strength of the left insulator in the system. As the ef- 
fective barrier strength changes, one sees negative cross- 
correlations turning positive for L < £. This indicates 
that a gate voltage can tune the entanglement proper- 
ties. More interesting is the case L = £, where noise 
cross-correlations turn completely positive enabled by the 
strong CAR signal. In the specular regime the noise is 
enhanced. This can be understood from Fig. [3] where 
the CAR signal in the specular regime is double than 
that of the normal case. The behavior depicted in Fig. [5] 
is of significance for the experimental detection of en- 
tanglement in solid state systems. It shows that a gate 
voltage can change the sign of noise cross-correlations 
unlike that predicted for normal metal counterparts. It 
is worth mentioning that for L » ( the magnitude of 
the noise cross-correlations are very much reduced (not 
plotted) but one can also see completely positive noise 
cross-correlations. 



IV. CONCLUSIONS 



in case of enhanced positive noise cross-correlations, in- 
dicating bunching. We point out the novel phenomena 
of specular crossed Andreev reflection, which enhances 
noise cross-correlations. The settings envisaged in this 
work are experimentally accessible. A typical s-wave su- 
perconductor like Aluminium has a coherence length of 
£ = 1600rtm. Since the proximity effect induces super- 
conducting correlations in graphene, an Aluminium su- 
perconductor on top of the graphene layer would give rise 
to a similar correlation length. This separation would not 
be a challenge since crossed Andreev reflection measure- 
ments are carried out routinely at less than these lengths. 
Further, the superconducting gap in Aluminium is lmeV, 
while the typical Fermi energy in normal doped graphene 
is around 80meV. In our study we have considered for 
certain situations Ep = 100 A, i.e., Ep ^S> A, which cor- 
responds to undoped graphene, while Ep <C A can be 
tuned via doping graphene. These values are realistic 
and thus obviate any reasons for scepticism. Employing 
these entangled states for quantum information process- 
ing will increase the allure of graphene. 



Recent CAR experiments^ are the next generation 
in detecting the splitting of Cooper pairs into differ- 
ent leads, thus probing entanglement in the context of 
nanophysics. In this work we provide the results of noise 
cross-correlation spectra as a function of gate voltage for 
a NISIN graphene junction. The Fano factor (not pre- 
sented here) is also on predictable lines and shows a spike 
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